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Abstract 

We calculate the decay rates <fi — > 77r°7r° and jit !] in very good agreement with the recent accurate 
experimental data. We also point out the necessity of a <fi-fK°K a contact vertex beyond the pure 
K + K~ loop model. The decay widths <f> — > 7/o(980) and <j> — ► 7ao(980) are also calculated taking into 
account the finite widths of the scalar resonances /o(980) and do (980). The latter are shown to be 
essential in order to obtain meaningful results. The resulting decay rates to 7/ (980) and 7a (980) are 
in good agreement with the experimental ones without invoking isospin breaking in the couplings of 
the /o(980) and ao(980) resonances to the K + K~ and K°K° channels, at odds with recent proposals. 
The derived formula for calculating these <j) radiative decay widths can be also applied in their own 
experimental analyses in order to obtain more precise results. 
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1 Introduction 



In this article we study the 4> radiative decays to 7 and two neutral pseudoscalars and the related rates 
to 7/o(980) and 700 (980). These decays were first observed experimentally by the CMD-2 [pi and SND 
[§, ||, U] collaborations at Novosibirsk by measuring the invariant mass distributions dT((f> —* 77r°7r°) /dm^ 
and dT((p — > jir°r]) / dm^ , respectively, where we indicate by ^ro the invariant mass of a pair of neutral 
pseudoscalars M° and N°. On the other hand, high statistics results have been recently reported by the 
KLOE collaboration at DA<£NE in refs.[j|, ^] for the rates <j) — > 77r°7r° and 77r r/, respectively. These decays 
offer an additional source of experimental information on the non-trivial and so important low lying scalar 
mesons with vacuum quantum numbers ++ . Our theoretical study generalizes the one undertaken in 
rcf.Q about the <j) ~^ 'yK°K° decay in order to reproduce the precise data of refs.|5|, |] on <p — > 77r°7r° and 
77r°?7 and connects the previous processes with chiral symmetry and unitarity. Ref.[7]] was also followed 
closely by ref.[]|] where the aforementioned radiative decays to the final states 77r°7r° and jtt°t] were also 
considered. For recent accounts on the application of similar techniques to the scalar meson- meson and 
meson-baryon dynamics we refer the reader to refs.|], |T(^, [TT|| w here many references are presented and 
discussed. Other studies of 4> -> jM°N° decays are refs~|2[||, ||, H, 0, 0, H, H, U H|. 

In order to interpret the experimental results from refs.|| |J on the 4> decays to 7/o(980) and 7ao(980), 
ref . |23| has pointed out recently the necessity of considering important isospin breaking corrections in the 
couplings of the /o(980) and ao(980) resonances to the K + K~ and K°K° channels. That is, 
and |fl^o/fol> where R represents either the /o(980) or ao(980) resonances, are no longer equal as required 
by isospin symmetry but can be actually quite different according to ref.|2^]. Special emphasis is given to 
the rather large decay width <j) — > 7/0 (980) together with the deviation from one, by around a factor four, 
of the ratio Br(<f> —* jfo)/Br(<f> — > 700) as experimentally established in refs.|l], ||], ]f| |||. Nevertheless, 
in ref.[[24|] it is clearly shown that possible isospin breaking effects in the quotient T((p — > 77r°7r°))/r(0 — > 
77r°?7) cancel to a large extend. As a result, this ratio is sensitive mostly to the isospin eigenstates 
ao(I = 1) and fo(I = 0). Similar conclusions are obtained within the vector meson dominance model of 



ref.|25|] as well. Other references treating the /o(980) and ao(980) mixing are |2(| ^7|, [2g, 29]. 

We show that one does not need to abandon isospin symmetry in the couplings of the scalar resonances 
/o(980) and ao(980) to the KK channels so as to account for the experimental results on the <fi — > jR 
decay widths. We emphasize the fact that the threshold of jR is so close to the mass of the 0(1020) 
resonance, that meaningful results for the width T((j) — > "fR) can only result once the finite widths of the 
/o(980) and ao(980) resonances are taken into account. It is important to stress that the decay width 
T((fr — > jR) depends cubically on the small photon momenta and hence small changes in the nominal 
masses of the /o(980) or ao(980) resonances, e.g. as compared to their widths or to the difference between 
the real parts of the pole positions and the peaks of the scattering amplitudes, imply dramatic variations 
on the resulting (f> radiative decay widths to jR if the standard two body decay formula were used. 
In sec.^ we calculate the </> decays to 77^77, 77r 7r° and 7^°^° where we connect with and extend the 
results already given in refs.0, |||. We also establish a new contribution beyond the pure K + K~ loop 
model of ref.(l^] due to a contact (f)jK°K vertex. In sec.^, making use of the results of the previous 
section, we take into account the finite width effects of the /o(980) and ao(980) resonances on the decays 
<p — > jR and derive the corresponding formulae. We also consider the connection of the obtained formula 
for T((j) — > 7i?) to the invariant mass distributions dT(<p — > 7 M°N )/ dm m°n°- Finally, our results are 
discussed in sec.|] and the conclusions are given in sec|| 
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2 T((f) -> 7 M°iV ) rates 

We are interested in studying the <f> — > 7M°iV decays. The final state interactions (FSI) of the two 
pseudoscalar mesons are assumed to be dominated by their strong S-wave amplitudes which contain the 
/o(980) and ao(980) scalar resonances. Then gauge invariance requires, as noted in refs. ]l2], |37], [l5| with 
respect to the related decays 4> ~~ > 7-R as well, that the transition amplitude ^(1020) — ► jM N must 
have the form (?], ||]: 

M [HP) - 7(k)M°N°(Q)} = H{Q 2 ) [g^{p • k) - p a k?~\ e* (7)^)0 . (2.1) 

The kinematics is indicated in fig.[l] and the four-momentum Q corresponds to the total one of the pair 
of pseudoscalars M°N°. 

As in refs. [12, 16, fjj we consider the K + K~ meson loop contribution to the decay <fi — ► jM°N°, fig.|l], 



since the K + K~ system strongly couples to both the cj> and R resonances. In addition, since the K°K° 
couples with the same strength as the K + K~ state both to the scalar and 4> resonances, we also allow 
for a neutral kaon loop contribution, as indicated in fig|l|a, where the vertex on the left of the figure is 
discussed in length below. 

To simplify the discussions, we will consider the isospin channels separately, since they do not mix 
by final state interactions in the isospin limit, so that H 1 (Q 2 ) is the invariant function in the same 
isospin channel as the corresponding scalar resonance R from the kaon loops of fig.|l[ We now calculate 
the function TL^Q 2 ) defined such that H 1 (Q 2 ) = Ti 1 (Q 2 )tK+K-->M°N ' with ijc+K--*M N° * ne i sos P iri 
amplitude — ixj^^MN ^m^n /^' wnere the KK and MN states are in the isospin channel I of the 
scalar resonance R, namely, I = for the /o(980) and I = 1 for the ao(980), and —l/y/2 and C 1 ^^ are 
Clebsch-Gordan coefficients. For the 77^77 final state only the 1 = 1 channel takes place and the 1 = 
one is the only contribution to 77r°7r°. For the jK°K° final state both contributions sum up. 

In ref.|?]] we took as interacting Lagrangian, 

C = 2e g<t) A^^K + K- - i (eA^ + g<t> ^)(K' d»K+ - d»K-R + ) , (2.2) 
upon making the (j) an d K + K~ interactions gauge invariant. 

7K + , K + K + 

,->«-fc^ y 

K , K° N ° K K S K - 

a ) b) c) d ) 




Figure 1: Set of diagrams for the calculation of <fi — > -fM°N° through kaon loops. 



We now consider the more refined approach of ref . [ 30 1 , where the couplings of the octet of pseudoscalars 
to the octet of vector resonances are given at order p 2 by taking into account chiral symmetry in a chiral 



3 



power expansion. The couplings of the singlet vector resonance u\ cancel at this order exactly. Making 
use of ideal mixing we can write the (ft resonance as: 

2 1 

(ft = —j=wa + -j=W. , (2-3) 

and thus the couplings of the (ft resonance from ref. fjOfl are the ones of the u>% resonance times — 2/\/6. 
The tree level matrix element for the contact vertex (ftjK + K~ , corresponding to the left vertex of fig.|l]a, 
contacts calculated from the formalism of ref. ||30| is: 

M((ft - 7 K+K-) contact = - ^ e ^f V e( 7 ) ■ e((ft) - j^L - Gy) Pa e((ft)p (k^f - fc* c ( 7 )«) , 

' _ ' . (2 - 4) 

where / ~ 93 MeV is the weak pion decay constant and Gy measures the strength of the transition 

(ft — ► K + K~ . The experimental width T(<ft — > K + K~) is reproduced with Gy = 55 MeV. Let us note 

that the second structure in the equation above is gauge invariant by itself. On the other hand, vector 

meson dominance requires Fy = 2Gy [^0| and in this limit the coefficient in front of the second term of 

eq.([2.4|) is zero. At the chiral order considered in ref.|3(]], there is no direct coupling (ft'-fK K although 

at higher orders this is no longer the case and can mimic e.g. (ft^K^K® vertices generated through the 

exchange of vector resonances, see e.g. [jig]. 

We consider first the K + K~ loops of figjljthat are originated from the surviving local term of eq.(|2^ 

fig.|l]a, plus the Bremsstrahlung ones, diagrams [l|b and c, and fig.|l]d. We follow the treatment of ref.| 

where it is shown that the T-matrix element t^ +K _^ K+K _ of ref.|3l| factorizes on-shell in these diagrams. 

The reason is, as proved in ref.0], that the off-shell parts of the T-matrices of ref. do not contribute 

to the coefficient of the term in eq.( |2.1[ ) proportional to p a k@, which is the function H(Q 2 ) as required 

by gauge invariance.p^ 

Motivated by the self- invariant structure of the second term of eq.( |2,4| ), we consider for each isospin 
channel a general contact (ft^KK interaction with the same structure but parameterized by a coupling 

Ci, 

ViicfryKK) = -^L Pa e(cft) p (*« £ ( 7 )/» - hPeiriT) ■ (2-5) 

The dependence of Q on the isospin channel implies a possible non-vanishing (ft^/K K local term since 
the Clebsch-Gordan coefficient for K°K° changes sign from / = to I = 1 while that of K + K~ is the 
same for both isospins. Similar structures for 77r°7r° and 77r°?7 are not considered since they are expected 
to be suppressed by the OZI rule. 

We are now in position to evaluate the loop of fig.|| with the contact term Vi(cjr)KK) of eq. (|2.5| ), in 



the vertex on the left of the diagram. The strong amplitude tK+K--*M°N° °^ re ^-@l (^k°k ^m°n ^ s j us * 
proportional to the former) factorizes on-shell since the off-shell part, as in the case of the pure strong 
interacting problem treated in detail in ref.[[n]], just renormalizes the coupling on the left part of the loop, 
namely, the Ci factor of eq.([2.5|). This is due to the fact that it only gives rise to tadpole like diagrams. 
Since the (j factor is a free one this renormalization process is not relevant for our present purposes and 
is just reabsorbed in the final value of £/. Since both vertices on the diagram of fig.§ factorize with their 
physical renormalized values, the one on the left does not involve any KK momentum, see eq. (|2.5| ), we 



#2 It is shown in ref. [[Io| that the contributions from crossed channels exchanges of pairs of pseudoscalars and resonances 
tend to cancel in the S-waves with / = 0, 1 and 1/2 below roughly 1 GeV. 
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Figure 2: Diagram for <fi — > jM°N° with the contribution of eq.(2.5), which is gauge invariant by itself, in the 
vertex on the left. This contribution was considered in ref.Q to some extend but not in r cfe|§0 [[!. Note 



that in the vertex on the right the full amplitude t% 



K+K-->M°N° 



is considered. 



are then left with the logarithmic divergent loop integral, 

A 1 



(2vr) 4 (g2 _ m 2 K+ + i0 +) ((Q _ g)3 _ m ^ + iQ +) 



(2.6) 



In ref.||] this integral was taken the same as the one fixed in ref.[^] for the strong interactions. Never- 
theless, one should point out that unitarity and analyticity, with the latter restricted to the presence of 
the unitarity or right hand cut, allow the presence of an extra subtraction constant, 5G 1 , in the integral 



of eq. fl2.6| ) as compared to G K ^_ of ref.[31]. We show now this point. 

Since the vertex of eq. (|2.5| ) is gauge invariant by itself, we consider its contributions separately. Thus, 
it is characterized by the corresponding invariant function H 1 , a contribution to the total one H . Let us 
denote by the invariant function corresponding to the <fi — > decay where a indicates any channel 
that can couple strongly to the KK channel and with the same well defined isospin. We treat these 
invariant functions as form factors depending only on Q 2 since this is the dependence obtained from 
unitarity loops as those of figj2| that we want to analyze]^] Then unitarity requires: 



Imag 



Hi] = Y.^)* pp ^— °(q 2 - w D > ( 2 - 7 ) 



where p@ is the phase space factor of channel 0, pp = qp/{&Tt\J Q 2 ). In ref.[32, 33] it is shown that when 
the T-matrix appearing on the right side of the previous equation has no crossed cuts due to crossed 
channels, as the one of refs.|3l], [34|, H^, collected in the vector column H 1 , can be represented as: 

H 1 = [1 + K\Q 2 ) ■ G{Q 2 )] - 1 RdQ 2 ) , (2.8) 

where it is important to remark that Rj(Q 2 ) is a vector of functions without any cut. In addition, 
T 7 (Q 2 ) = [1 + K 1 {Q 2 )G{Q 2 )]~ l ■ K^Q 2 ) with T 7 (Q 2 ) the strong T -marix for the S-wave meson-meson 
scattering with isospin /, or 1. In ref.plj K 1 (Q 2 ) corresponds to the lowest order CHPT meson-meson 



amplitudes [g^, 36] while in ref.[34] together with this contribution one has the local terms calculated from 
0(p 4 ) CHPT, see refs.pl|, |3~i| ] for explicit formulae. In addition G(Q 2 ) is a diagonal matrix of unitarity 
loops such that for 1 = 0, Gn{Q 2 ) is the loop analogous to eq.(^l]) but for 7T7r(a = 1) and G22{Q 2 ) 
corresponds to the KK(a = 2) channel. Analogously for 1 = 1, G\i{Q 2 ) corresponds to -Krj{a = 1) and 
G 22 {Q 2 ) to KK{a = 2) and is equal to G 2 2{Q 2 ) with 1 = 0. 

On the other hand, because of the aforementioned factorization of the T-matrix of ref.[]3l|], H^ can be 
expressed as: 

H^ = \i-T I g I (Q 2 )\ I , (2.9) 

#3 Note that the vertex on the right is the full strong T-matrix with its own unitarity cuts. 
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with Aj the vector column of direct couplings to (jryKK after removing the tensor structures like in the 
definition of H 1 , that is (A/) T oc (0, (i). On the other hand, Qli{Q 2 ) can differ at most from Gu(Q 2 ) in a 
subtraction constant which we express by writing Q 1 = G + 5G 1 . Notice that because (A/)i = only the 
matrix element of 5G 1 for the KK channel is relevant, the only one that we will keep in the following. 
Now, taking into account the explicit expression of T 1 in terms of G and K , one has: 

H 1 = (l- [l + K^Y 1 {^G + K^G 1 }^ A/= [l + i^Gj^A/- [l + K 1 G]' 1 K 1 SG 1 Xi 

= [l + K^Y 1 (l-KHG 1 )^ , (2.10) 

and hence Rj = (1 — K 1 dG^Xj. Let us note that Ri has no cuts, as it should, since neither K 1 , 5G or 
Xj have any cut. 

As a result of this digression, we will keep 5G 1 as free parameters in addition to £j. This amounts to 
four free parameters, two for each isospin channel, although in the end we will show that two of them 
will turn out to be fixed. Although eq. fl2,9| ) is derived by taking into account the off-shell form of the 
T-matrices of ref . |3l| , we will also use it for the more involved T-matrices of refs.[^] since the equivalent 
form given in eq. Q2.10|) is equally valid for all of them. That is, the column matrix Rj(Q 2 ) is fixed to be 
(1 - K^G^Xt, eq. (|2.10D , for both the T-matr ices of ref.|3l| and ref. ]j34j| . While for the former strong 
amplitudes this expression for Rj is a strict derivation due to its specific off-shell parts that enter in the 
unitarity loops, for the latter ones Ri is fixed by analogy. 

Summing the K + K~ loop contribution of figs.|l]a, b, c and d plus the KK loop of figj2|, H 1 , we can 
then write: 

= l f:^ V I(a,b) - ^-QiAQ 2 )]^-^ , (2.11) 



where for the 77r°7r° {i^rj) final state only 7 = 0(1) contributes while for ^K^K® one has to sum both 
1 = and 1 contributions together. From this expression it follows that Ti 1 (Q 2 ) = H 1 (Q 2 )/t^- +K _^ M0N0 
is given by: 



»'<« 2 ' ^ s^f>»> " ■ (2 ' 12> 



Including the three-body phase space and taking into account eq.( |2.11| ), the width T(<j) — > ^M®N°) for 
M°N° = 7r°7r° or ir°r] can be expressed as: 



r(0^ 7 M o Af°)= dy/Q' 



)2 



q|k| 3 |pM| 




M <t> G v T(n h \ y^a 1 (o 2 u t R 



(2-13) 

where |pm| is the momentum of the M°N° system in their center of mass frame and a is the electro- 
magnetic fine structure constant. Possible symmetric identity factors 1/2 are omitted in eq.( |2.13| ) since 
the transition matrix element | t^ +K -_^ M o NO | 2 is taken from refs.|Jl], 34] and such factors are already 



included in its normalization. In the previous expression the sum is restricted over those isospin channels 
included in M°N°. For the ^K°K final state case, M°N° = K°K°, one has to add to the term between 
bars in eq.( |2.13| ) the contribution +(Co — Cl)/v / 2-^<4 from the direct local terms of eq,( |2.5| ) with 1 = 
and 1. 
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3 r(0 — ► jR) decay widths 



As in the previous section, see also refs.|l2|, [37], |15|, |7|, gauge invariance requires that the transition 



amplitude 0(1020) — ► jR must have the form: 

M [HP) - 7(k)R(Q)] = H R (Q 2 ) [g^(p ■ k) - p a k?} e A (7)a^)/3 • (3.1) 

where now Q is the total four-momentum of the resonance R. The function H R (Q 2 ) can be obtained from 
the expression given above for (Q 2 ^ replacing ^^■+/^-_ > jy,fo^yo hy the coupling of the K~^~ K channel to 
the R resonance, g R +K _. That is, we are extracting the resonance pole contribution from tx+K-^M°N°- 
Although not explicitly shown g R +K _ can depend on the energy, Q 2 , flowing to the resonance R which 
is not fixed since its mass is not well defined due to the finite widths of R. Indeed, the latter gives rise to 
important effects due to the proximity of the mass of the 0(1020) resonance to the nominal ones of the 
/o(980) and ao(980). Let us consider this in detail. 

Since the coupling g R + K ~ appears as a factor in H R (Q 2 ), then we have that H R (Q 2 ) = Ti 1 (Q 2 ) g R :+K -- 
It is also convenient to define the quantity M' so that M(<p — > jR) = M.' g R +K _. Then, the width 
T(4> 7-R) can be expressed as: 



v — 1 A — 1 

Of course, because of three-momentum conservation, Q = — k in the center of mass frame where the 4> 
is at rest. We have included energy distributions for the scalar resonances /r(Q°) in order to take care 
of their finite widths. This is not considered in the literature for the calculation T((j) — > 7I?) since well 
defined values for the masses of the scalar resonances /o(980) and ao(980) are assumed and the standard 
two body decay formula, eq.( 3.10| ), is then applied, see e.g. (l^, ^5|. We should note that this is a 



very unstable procedure since the decay width scales as |k| 3 , where |k| is the photon three-momentum: 

... M l~Q 



2Ma. 



(3.3) 



Allowing a change of rriR of just 10 MeV from the central value mR = 980 MeV, the quantity |k| changes 
by a factor of 1.7 and |k| 3 by a factor of five. The standard two body decay formula can be obtained 



from the more general eq.(pl) by taking f R {Q°) = 5(Q° - JQ 2 + ml). 



In order to fix Jr(Q°) we invoke unitarity which implies: 

2Im ag [t R +K ^ K+K .) = {K+K-\T^T\K+K-) . (3.4) 

Now, assuming that unitarity is saturated by the exchange of the corresponding scalar resonance R, it 
follows that, 

2l™ g [t R K+K ^ K+K A = f j0^f R {q"){2nf5\Q- q )\ g R +K .\ 2 

= ^o/r(<5°) \g R + K -\ 2 , 



(3.5) 



and then, 



^ofdQ°) \9k + k- I 2 = l^g[t^ +K -^ K+K - (Q 2 )} ■ (3.6) 
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The previous equation gives our final form for fn(Q ). We show in fig.|| the combination 7rfn(Q )/2Q = 
lmag[t^ +K _^ K+K _ (Q 2 )] for the /o(980), left panel, and ao(980), right panel. The dashed lines corre- 
spond to the matrix elements t^ +K _^ K+K _(Q 2 ) given in ref.|3l]], where the only one free parameter is 
fixed in this reference by a successful study of the strong tttt, KK and ttt] S-wave meson- meson inter- 
actions. These T-matrices have been used by now in a whole series of studies about the scalar sector 
stressing the role of unitarity and chiral symmetry in order to understand the scalar dynamics. These 
studies comprise strong interactions [31], fusion of two photons to 7r°7r°, 7r + 7r~, K + K~, K°K° and 7r°r/ 



[38 1, (j> radiative decays 0, J/^ decays [32] and pp collisions M. In the same figure we also show by the 
solid lines, Imag[t^ +K _^ K+K _ (Q 2 )] from the T-matrices of ref.[34] with the use of the Inverse Amplitude 
Method p9[ . We also compare our choice for 7t/r(Q°)/2(5 from the T-matrices of refs.[31, 
one that results by taking, 



Ml, with the 



V R {Q 2 



-Imag 



Q 1 



m. 



v r (q 2 ) j 

> - Ren(m|) + n(Q 2 ) 



(3.7) 



from ref.[12], since this resonant form was followed in the experimental references jl], |2|, ||, |], || || to fit 
their data. The energy distribution resulting from eq.(|3,7j) is shown in fig.|3] by the thin dotted lines. The 
free parameters in eqs.([T7]), two masses and four couplings, are fixed to the central values given by the 
best fits of refs.|3|, Q to dT(4> — > 77r°7r°)/dm 7r7r and dT((j) — > 'yTr r))/dm 7rrj , respectively. The Ur(Q 2 ) is 
the sum of two-meson self energies, tttt and KK for 1 = and tttj and KK for I = 1 as in refs.[||, 
Explicit formulae for Il(Q 2 ) can be found in ref.[12|. It is interesting to remark the presence of a long 
tail in each of the energy distributions fn(Q°) towards low invariant masses. This is a result of keeping 
the real parts of the unitarity two-meson loops |31, [?4|] or the real parts in the two-meson self energies in 
V R (Q 2 ), eq.Q 




(MeV) 




(MeV) 



Figure 3: nf R (Q°)/2Q = Jm&g[t% +K _^ K+K _} for the /o^80), left panel, and a (980), right panel. The solid 
and dashed lines are calcul ated from the T-matrices of refs.|34|, [Uj, respectively, and the thin dotted ones from the 
resonant form given in eq.(3.7) following ref. |]l2| and the experimental analyses |J. The values of the masses 
and couplings, necessary to apply eq.(|3.7|), are determined by the best fits of the latter references. 



Then, T( 



jR), eq.(|3.2|), can be written from eqs^ST) and ( |3.6| ) as: 



2\2 



(3. 
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Of course, one should keep in mind 



where we have used the shorter notation instead of ijc+K~-^K+K 

that Q 2 and |k| are not independent variables but related by energy-momentum conservation, eq.^ 
It is also apparent from the previous equation the cubic dependence on the photon three-momentum 
From eq.(|3.8|) we can also define the differential decay widths, 



dT{4) -» ryR) _ M$ dT((j) -► jR) 



d\k\ 



J^Llmag [tg(Q 2 )] \n\Q 2 )\ 2 (M 2 - Q 



2\2 



(3.9) 



Note that Imagt^ extends from the two pion threshold for the /o(980) case and from the nr] one for the 
ao(980) resonance. 

As stated above, the standard two body decay formula for (j) — > jR, with a well defined mass m R for 
R, can be obtained from eq^^) by performing the substitution = 6(Q° — wr). In this case the 

decay width is given by: 

T{ ( f ) ^ 1 R) mRflxed = ^\H I {rn 2 R )gK +K -\ 2 . (3.10) 

Following the previous lines, we can also derive an expression for the decay widths T(cj) — > jM°N°) 
by considering explicitly that the strong amplitudes tjc+K--*M fi N° are dominated by the exchange of the 
corresponding R resonances. One should keep in mind that the results given in sec.§ do not make such 
assumption and are derived directly in terms of the corresponding strong amplitudes t^ +K _^ M0N0 . 

Let us denote by qi and q2 the three-momenta of the M° and N° particles and by wi and w 2 their 
corresponding energies. Then, since we are considering the decay <fi — > jM°N° to be mediated by the 
resonance R, if follows that, 



at0\ 



7 M u iV 



3 2 



dk 



(2vr) 3 2|k| 



Ml 



V R (Q 2 



2VQ 2 



2^/Q2j (2tt) 3 2w : 



dpx 2ir5(Q° - wi - u> 2 ) , R |2 
— \9m°n°\ 



2w 2 



(3.11) 



where T> R {Q 2 ) is any possible inverse propagator for the resonance R. Let us note that the term between 
curly brackets just corresponds to the standard formula for the decay width of a resonance of mass y/ Q 2 
to the M°N° state with a coupling g^p N o- In addition, we can rewrite, 



M 


2 


M'g* +K „(Q 2 ) 


Vr(Q 2 ) 




V R {Q 2 ) 



\M'\ 2 |4(Q 2 



(3.12) 



On the other hand, within our present assumption that t^iQ 2 ) is dominated by the exchange of the 
scalar resonance R, unitarity implies: 



Imag[4(Q 2 )] = ^^Y ^f^ ^(Q 2 ) = |4(Q 2 )| 2 El^! 2 ^ 



wl 



+R 1 2 
c 22l 



R,tot 



4TTy/Q 2 



SttJQ 2 



(3.13) 



Where q a is the center of mass three-momentum of the two-body channel a, W a is the corresponding 
threshold energy and T R t t{Q 2 ) is the total energy of a scalar resonance of mass \J Q 2 with the same 
couplings, indicated by g^, as the R resonance. 
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As a result of eqs.( |3.12| ), ( |3.13 ) we can then write instead of eq.( 3.11 ): 



atO\ 



jM N 



3 2 

-YY — 

v=\ A=l 



^l^flmag[4(Q 2 )]^^ 



(27T)3|k| 1 ' ' ° L ^ ,J IWQ2 

2 ^l^(Q 2 )| 2 Imag [^(Q 2 )] (M 2 - Q 2 ) 2 Br^ l0NO (Q 2 ) , 



(3.14) 



where r RA/ /Ojvo(Q 2 ) is the decay width of a scalar resonance with mass \J Q 2 and the same coupling 
(Q 2 ) to the M°iV channel as the R resonance at that energy. This parameterization can be seen 



9m°n° 

as an alternative to that of ref. [12| and is specially suited for those approaches that directly work with 
the strong amplitudes, without considering specific forms for the propagators of the scalar resonances 



which are explicitly needed in the legitimate parameterization of ref. 1 12]. 

Consistently with eq.( 3.13 ), the branching ratio Br^, l0N0 
calculated as: 

I R 1 2 

M fr>2\ — \9m°N°\ Pm°n° 



(Q ) of the strong decay R(Q Z 



M°N° 



is 



Br 



M°N° 



(Q 



(3.15) 



where in addition one should include, when necessary, the corresponding 1/2 factors for the decays to 
symmetric two particle states. For the sake of completeness we show in table |l| the pole positions, \/~sr 
and couplings, g^, for the /o(980) and ao(980) resonances from the T-matrices of ref. pl| . ^ The latter 
are defined by: 

SR)t^(s)\ , (3.16) 



\9a9p\ 



lim |(s 

s^s R 



where sr is the pole position of the R resonance in the unphysical sheet where the center of mass three- 
momentum of the M°N° state has negative imaginary part, while that of the KK three- momentum is 
positive. It is also worth noticing that the couplings |<?5^-| are different for the /o(980) and ao(980) 
resonances like in ref. |l9| . 

Taking from table |l] the values of the couplings and pole positions one calculates from eq.( |3.15| ), 



Brl° w = 0.7 and Br^ = 0.63 in perfect agreement with the branching ratios given in ref.[31] where the 



finite widths of the /o(980) and ao(980) are taken into account for the strong decay rates as well. 



/ (980) 


a (980) 


V /Sfo = (992.6 - i 11.8) MeV 
\gt \ = 1.90 GeV 
\gfc R \ = 3.80 GeV 


^/s^ = (1009.2 - i 56.1) MeV 
Isfi* | = 3.54 GeV 
\g% R \ = 5.20 GeV 



Table 1: Poles positions, 
matrices of ref. [31]. 



y/sR, and couplings, g&, of the /o(980) and ao(980) resonances from the T- 



Comparing eq.( |3.14| ) with eq.( |3.8[) for the decay <j> — > jR one obtains the relation, 

djkj " djkj Br M°N0(Q ) > (3-17) 



* 4 Since the 1 = tt-k channel is completely symmetric under the exchange of the pions, the convention of including an 
ext ra fac tor l/y/2 in the definition of this state is followed in ref. blfl. Thus, the 7T7t coupling resulting by applying directly 
eq.(3.16) to the 1 = T-matrix of ref. has been multiplied by y2 in table [l| 
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which, together with eq.( |3.9| ), can be used to analyze in a pure phenomenological way the experimental 
data of (ft — > 7M°iV by parameterizing and fitting Im&g[tx+ K -^ K + K -], as e.g. in eq.(3.7). Because 
we are assuming that the decay (ft — > 7M°iV is dominated by the exchange of the corresponding R 
resonance, the previous relation tells us that the differential width (ft — > jM N is the same as the one 
to 7-R multiplied by the probability that the scalar resonance decays to M°N°. The latter is just the 
branching ratio since the resonance will certainly decay to any channel at asymptotic times. This natural 
result is indeed the way the experimental values on the (ft — ► 7/0 (980) and (ft — ► 700 (980) widths are 
obtained once the rates (ft — ► 7/0 (980) — > 77r°7r° and (ft — > 700 (980) — ► 77r°?7, respectively, are determined 
by fitting the experimental data Jl], |3|, |||, ||.^ Nevertheless, one should keep in mind that eq.( |3.17| ), 
when used with the simplified form of Br^ [0NQ given in eq. Q3.15 ), is an approximation since it does not 
take into account the finite widths of the /o(980) and ao(980) resonances in their strong decay rates. 

Let us stress that eq.( |2.13| ) is given in terms of the strong S-wave T-matrices of ref. [31, |34| without any 
reference to the exchange of an scalar resonance R. Indeed, in ref.pl[] the scalar resonances /o(980) and 



ao (980), together with the a meson, are generated dynamically as meson- meson resonances in terms of an 
interacting kernel given by the lowest order CHPT amplitudes 



36] without any explicit resonance field. 



Thus, it is an important consistency check whether the invariant mass distribution obtained from eq. (|2.13| ) 
agrees indeed with the one of eq.(|3.17|) when substituting in the latter dT((ft — > 7-R)/d|k| from eq.([3.9|). It 



is straightforward to demonstrate by relating Imag[i^- + ^- 
that eq.( pT3 ) and the combination of eqs.( |3.17 ) and |T 



and 



.K+K-i \"K+K— *m°N°(Q 2 )\ 2 via unitarity, 
yield the same dT((ft -> 7M°iV )/d|k| below 



the KK threshold. Above it, a specific energy dependent form of Br^ 0N0 should be employed. 



^. 40 

T 
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Figure 4: Invariant mass distributions dBr(<ft — > 77r°7r° )/c?m W7r 10 8 and dBr((f> — > jTr°rj)/dm 7rn 10 s from left to 
right, respectively. The thick lines correspond to a simultaneous fit of the parameters Cr and 5G 1 to the experimental 
points of H (empty circles) and (full circles), for the 77r°7r° final state, and to the data of Q (empty circles) 
and |6| (full circles) for the jir°r] one. The thin line corresponds to £/ = —Gy/VS and SG 1 = 0, the appropriate 
values to reproduce the results of ref. B . 



#5 In ref.Q an extra and large 0-7 contribution is included in addition to the fo (980)7 one with a destructive interference 
at low 7T7T invariant mass. However, this is a controversial result at odds with the conclusions of the previous experimental 
analysis refs.^, |^|, giving rise to an unconventionally large rate <f> — + 7/0 (980) larger than the ones given from any 
model. 
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4 Results and discussion 



In sec.§ we have derived the corresponding expressions for the <f> — ► jM^N decay by taking into account 
the final state interactions due to the strong S-wave meson- meson amplitudes. This is what we denote 
by the scalar contribution to such decay rates. Nevertheless, for the 77r°7r° final state the background 
4> — ► vr°/3 — > 77r°7r° is not negligible, its relative size depending strongly on the energy region, as we show 
below in fig]5|. In our final results we have included the interference term between the scalar and pir° 



contributions with our own scalar amplitudes and with the vector part calculated as in ref.[22|. For the 
appropriate formulae for the vector piece we refer to that reference. 

We now perform a simultaneous fit, including both the scalar and vector contributions, to the 4> — * 
77r°7T data of refs.||, || and those of <p — ► jir°r] of refs.Q, r|. The fit to the invariant mass distributions 
dT((j) — ► 77r°7r°)/dm 7rjr and dT(<j) — ► 77r c '77) / 'dm^ with the T-matrices of refs.|31, 34] is presented in figj^in 
the left and right panel, respectively. The solid lines corresponds to the final state interactions calculated 



from the scalar amplitudes of ref.[34] and the dashed ones to the amplitudes of ref.[31|. The fit reproduces 



fairly well the experimental invariant mass distributions and the resulting values of the parameters are: 

T-matrix Co = +164.12 MeV 5G° = 1.46/(4vr) 2 
ref.Q Ci = -165.87 MeV 5G 1 = 1.36/(4^) 2 

T-matrix Co = +124.99 MeV 5G° = 1.61/(4vr) 2 
ref.|l Cl = -132.26 MeV 5G 1 = 1.44/(4tt) 2 

Table 2: Values of the parameters from a simultaneous fit to dT(cj) — ► r yir 7r°)/dm. K7T ||, || and dT((j) — > 
jn°ri)/dm^ |,|]. 



From the previous table of values is clear an obvious symmetry in the results. We have obtained 
that 5G° ~ 5G 1 and Co — — Cl- The Cr parameters represent the direct coupling (p^fKK with isospin 
/, eq.( |2.5| ). Hence the direct coupling for the K + K~ channel as given by eq.([2.5|) and the values in 
table |2| is proportional to — (Cl + Co)/\/2 ~ and vanishes. On the contrary, the one for K°K° is 
(Ci — Co)/v2 ~ \/2Ci- Let us note as well that because 5G° ~ 5G 1 one has that Q l ~ Q°. Since the 
Clebsch-Gordan coefficients for K + K~ are the same both for 1 = and 1 its contribution to any process 
is just proportional to £?°Co + G 1 (i — and cancels, as one should expect since the direct coupling 
(j)'yK + K~ on the left vertex of the diagram of fig.|| vanishes as we have just seen. Then it is clear the 
consistency of having obtained within the isospin formalism that Q 1 = G° once Co = — Ci) otherwise there 
would have been a mismatch between the results obtained from the physical basis of states and that 
of isospin. For the K°K° channels, since the Clebsch-Gordan coefficients change of sign when passing 
from J = to J = 1, one has the non-vanishing result Q 1 ^ — <7°Co — 2<5 1 Ci- Thus our results suggest 
the existence of a <p"yK°K° local term of the same type as that of eq. (|2.5| ) but with Cr replaced by 
Ck°k° = v^Ci an d the absence of such terms for K + K. The resulting value for |C/| is very similar to Fy 
from p — ► e + e~, Fy = 154 MeV or from 4> — ► e + e~, Fy = 165 MeV. The value for SG 1 is of natural size 
since it is a number of order one over 167r 2 . It is also important to remark that in the reproduction of the 
lowest energy part of the dT(cf> — ► 'jiT Tr )/dm nn invariant mass distribution, namely below 600 MeV, the 
background pir° plays a significant role, while its negligible above 700 MeV in agreement with ref. [|l9| . 
From the integration of the invariant mass distributions we obtain the branching ratios: 

Br{(j) -► 77r°7r°) = 1.0910" 4 , 
5r(0^77r°r ? ) = 0.72 10 -4 . (4.1) 
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for the T- matrices of ref . [ 34 1 . The values of the previous branching ratios when using the strong ampli- 
tudes of ref.[pl| are almost the same: 



Br((j) -> 7vr 7r°) 



Brl 



77T°?7) 



1.09 10" 4 
0.73 10" 4 



(4.2) 



7 7r 7r°) 



(1.09 ±0.03 



stat 



For comparison, the precise experimental results from refs.|| |(| are Br(( 
0.05 sys ) 10 -4 and Br(<fi — ► 77r°r/) = (0.796 ± 0.07) 10 -4 , respectively, in excellent agreement with our 
results. 

In fig.|| we present a new fit to the same data as before but imposing the constraints 5G° = 5G 1 and 
Co = — Ci that have emerged in a consistent way from the fit described above. The values obtained are: 



T-matrix of ref. [31] Co 



T-matrix of ref. [34] Co 



-Ci = +180.83 MeV , SG = 8G X 
-Ci = +146.42 MeV , 5G = 5G X 



1.42/(4tt) 2 
1.54/(4vr) 2 



(4.3) 



The quality of the fit is similar to that of fig]||. In the same figure we also show by the thin dotted line 
the pn intermediate contribution which as told above turns out to be relevant only in the lowest energy 
range of the spectrum. Let us note as well that in the experimental analyses at least 6 free parameters 
[|], ||, ||] are needed to fit the experimental (f> — > jtt ?] and 77r°7r° invariant mass distributions. 



^. 40 

T 

o . 
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Figure 5: Invariant mass distributions dBr((f> — > 77r°7r°)/(im W7r 10 8 and dBr(cf> — > ^ir^rf) / dm^ 10 8 from left to 
right, respectively. The thick lines correspond to a simultaneous fit of the parameters £/ and 5G 1 to the experimental 
points of H (empty circles) and jj| (full circles), for the 77r°7r° final state, and to the data of Q (empty circles) 
and |(| (full circles) for the ^yn !] one. The constraints 5G 1 = SG° and Co = — Ci have been imposed. The thin 
dotted line corresponds to the pn background. 

Once the Q and 5G 1 are fixed, table || and eq. (|4.3| ), we show in fig|], from left to right, the distributions 
dBr(cj) — ► 7/0 (980) ) /dm 10 8 and dBr(<j) — > 700 (980)) /dm 10 8 , respectively, calculated from eq.( p.9| ) and 
the strong amplitudes of ref.|3l], 34 1. As in figs.|| the solid lines corresponds to ref.[]34| and the dashed 
ones to ref.plj. We have denoted by m = \/Q^ the invariant mass of the /o(980) and ao(980) resonances. 



It can be surprising that the tails of the invariant mass distributions extend well below the prominent 
peaks of the /o(980) and ao(980) resonances towards rather low invariant masses. This is due to the 
cubic dependence on the photon three-momentum |k| = (M 2 — Q 2 }/1M§ which largely enhances the 
low energy part of the invariant mass distributions. Indeed, if we fix |k| to the value corresponding to 
some nominal mass of the resonances, e.g. mn ~ 986 MeV, then one obtains very peaked distributions 
around the masses of the /o(980) and ao(980) without any tail towards low energies. Hence, it follows 
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that although the resonance structure is clear and very prominent the low energy components of the 
energy distributions /r(Q°), see fig.||, cannot be neglected because of the enhancement due to the cubic 
dependence on the photon three-momentum, fig|| 
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Figure 6: Invariant mass distributions dBr(4> -> 7 / (980) )/dy / Q 2 10 8 and dBr(cf) -> 7a (980))/<i v / Q 2 10 8 
from left to right, respectively. The solid and dashed lines are calculated from the T-matrices of the 
refs.pl, 31], respectively. 



We integrate now the invariant mass distributions dT{4> — ► jR) jd\J Q 2 from the corresponding thresh- 
olds up to Q 2 = Ml, eq.(p3^), and we obtain the values shown in table |3[ 





T-matrix ref. 31] 


T-matrix ref. []34j 


Br(cf> -> 7 /o(980)) 10 4 


3.19 


3.11 


Br {</>^ja (980) ) 10 4 


0.73 


0.73 


Br(0— 7 / (980)) 
Br(<j>^-ya (980)) 


4.37 


4.26 



Table 3: Br( 



7 /o(980)), Br( 



7«o 



(980) and the quotient of both from eq.( 2.12| ) and the energ; 



distribution fn(Q ), eq.(|3.6D, determined with the T-matrices of refs.[34, 31]. 



In this table the energy distributions fn{Q°), eq.( |3.6D , are determined from the T-matrices of refs.pH 
|34"| ]. For Q = —Gy/VS and SG 1 = 0, which corresponds to the results of ref.||, one has Br(4> — > 
7/o(980)) = 2.34 10~ 4 , Br(<f> -> 7a (980)) = 0.86 10" 4 and the quotient between both rates is then 2.72. 

The SND collaboration, refs.|3| and Q, reports the branching ratios Br{4> — * 7 /o(980)) = (3.5 ± 
- 3 -al) 10 " 4 and Br O 7ao(980)) = (0.88 ± 0.17) 10" 4 , respectively. The CMD-2 collaboration, 
ref.Q, reports Br{(j) -» 7/ (980)) = (2.90 ± 0.21 ± 1.54) 10~ 4 . Taking into account the latter value for 
Br(4> — > 7/0 (980)), as in ref.jll]] since it arises from a combined fit to the 77r°7r° and 77r + 7r~ data^ 



and the one of ref.g for Br((j) -> 7a (980)), one has Br{<p -> 7 / o (98O))/.Br(0 -> a (980)) ~ 3.3 ± 2.0. 
These numbers are based in a parameterization of the experimental event distributions of 7r°7r° and 
7r°?7 leaving the masses and couplings of the /o(980) and ao(980) as free parameters. Once these are 
determined by the fitting procedure, the branching ratios to jR are obtained by dividing the resulting 
Sr(07 — > jR — > ^yM°N°) by the corresponding Br^ l0N0 . There is some model dependence on the 

#6 In addition, the error in the result Br((f> -» 7/0 (980)) = (2.90 ± 0.21 ± 1.54) 10~ 4 of ref.Q is enlarged in this reference 
so that the rate cj> — > 7/0, determined by performing a narrow pole fit, is also compatible within errors. 
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"experimental" numbers because of the specific forms of the resonance propagators. Nevertheless, we see 
remarkable agreement between these experimental results and our calculations of table |||. In ref.|| the 
KLOE collaboration reports the result Br((ft — > 7/o(980)) = (4.47 ± 0.21) 10 -4 . This value is obtained 
by including a very specific destructive interference in the low energy region between the fo (980)7 an d 
a o"7 contribution, both of them very large at such energies. The Br{(f> — > 700 (980)) is much clearer 
and very well established due to the absence of any significant background contributions to <j> —* jir°r] 
as established both theoretical and experimentally, refs.|L9|, [l7], 22, [|, g]. The KLOE collaboration Q 
determines the number Br((j> — > 700 (980)) = (0.74 ± 0.07) 10 -4 from the very precise measured rate 
Br((f> — > ■j7r 7]) = (0.796 ± 0.07) 10 -4 . Our calculations for both rates presented in eq. (|4.l| ) and in table ||| 
are in perfect agreement with these determinations since we described very accurately the corresponding 
invariant mass distribution as noted in figs.§ and |j| 





CMD-2-fit 1 U 

Br(cf> -> 7/0) 10 4 


CMD-2-fit 2 y 
Br{<j> -» 7/0) 10 4 


SND-g 
Br(4> -> 7/0) 10 4 


SND-|] 

Br((j) -> 7 a ) 10 4 


reported 


2.90 ±0.21 ±0.65 


3.05 ±0.25 ±0.72 


4.6 ±0.3^ 


0.88 ±0.17 


eq.(|2.12), Cl = 


3.21 


3.51 


4.8 


0.96 



Table 4: Br{4> — > 7/0 (980)) and Br((j) — > 7ao(980)) branching ratios where the energy distribution fn(Q°) 
is fixed from eq.( |3.7[ ) and then eq.( |2.12| ) is used. The results reported in the experimental references ||, Q 
are also given. In the second row of numbers, the function Ti 1 (Q 2 ), eq.( 2.12| ) is used with (j = 0, the 
choice that corresponds to refs.O, 0]. 



We now apply eq.( 2.12j ) to the energy distributions /r(Q°) coming from eq.(^), where the free 
parameters, the masses of the /o(980), ao(980) and four couplings, take the values of several fits of 
refs.[|], ||, |H used to obtain their experimental numbers of the rates 4> ~^ jR. The calculation is performed 
with the function H 1 (Q 2 ) calculated as in refs.(l|, |, |, |, that is, with (i = in eq. pTl2D . This 
corresponds to the second row of numerical results of table |||. In this way, the previous decay widths 
are calculated without the shortcut to use a fixed value for the branching ratio of the resonance R to 
the lightest decay channel in all the energy interval up to Ma, as in the present experimental analyses 

|, |, l-il- 
In refTg the quoted Br(<j> -> 7/o(980)) = (3.5 ± 0.3±J;|) 10" 4 comes by multiplying by three the 

measured Br(4> — > 77r°7r°) with Br^° (Q 2 ) = 1/3, see eq.( 3.17 ). In table |^ the quoted value from ref.Q 
is the one obtained with a fit to <fi — > 77r°7r° invariant mass distribution including as well a background 
from pir°, which amounts at most to 15%. The error is correspondingly enlarged so as to make both 
results compatible. For the ao(980) case the only number that is reported, both in table |] and ref.Q, is 
the same as the measured Br((f> — ► jtt rf) assuming Br^(Q 2 ) = 1. The fit 1 for the CMD-2 collaboration 
includes both the 77r°7r° and 77r + 7r~ final states. The CMD-2 second fit includes only the 77r°7r° final 
state. We see a general agreement between the reported and calculated numbers from eqs. fl3.7l ), ( Ijj ) 
and fl2l|), given in the first and second rows of numerical results, respectively, although the latter tend 
to be somewhat larger, particularly for the reported values of the CMD-2 collaboration. 

It has been recently claimed in ref.^J that in order to interpret the experimental results of refs.||, ||| 
for the rates <f) — > one needs to include sizeable isospin violating effects in the couplings of the 
/o(980) and ao(980) resonances to the K + K~ and K°K° channels, so that 7^ |9^o/fol an d the 

difference is claimed to be as large as a 30%. Particular emphasis is given to the necessity to deviate 
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from isospin symmetry in order to understand the quotient between the branching ratios of the (ft to 
7/0 (980) and 700 (980) that, although with a large experimental uncertainty, as shown above, has a 
central value of around three instead of one. The much more precise results of refs.0, 0], when taking 
Br(<p -> 7/ (980)) = 3Br((ft -> 77r vr ) and Br(<p -> 7a (980)) = Br(0 -> 7^%), imply a value 4.1 ± 0.2 
for this ratio, very close to our results of table ||[ Indeed in our approach the calculated Br ((ft — ► 7/o(980) 
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from eq.(2.12) and the one obtained by multiplying by three the Br(<ft — > 77T 7T ) differ in less than a 
Our study clearly shows that one can achieve a good agreement with the experimental data without any 
deviation from isospin symmetry in the couplings of the /q(980) and ao(980) resonances to the K + K~ 
and K°K° channels, although a contact interaction term (ft^K°K°, beyond the pure K + K~ loop model 
of ref.[12|, has to be included as described above. It is also clear that one should abandon in the study 
of the (ft — > 7-R decays the standard two body decay formula, eq.( 3.1C| ), used in refs.[23, 25], due to the 
proximity of the threshold of the final state to and the cubic dependence on |k|. As a result, the 
effects of the finite widths of the scalar resonances, and their associated energy distributions /ij(Q°), must 
be included from the very beginning. For instance, had we used the values for the /o(980) and ao(980) 
masses and the g^ +K ^ couplings given in table [l] we would have obtained Br(<ft — > 7/0) = 1.17 10~ 4 and 
Br((ft — > 700) = 0.58 10~ 5 . The latter value is so small due to the rather high pole mass of the ao(980). 
Let us stress that this pole mass is clearly different to the value where the S-wave 1=1 T-matrix elements 
peak, around 986 MeV, which indeed changes from one matrix element to the other. Hence, it is rather 
artificial to decide which is the value of to be used for the calculation of the rate (ft — > jR in an 
extraordinary sensitive two-body standard decay formula to the chosen ran value. 

It is worth mentioning that we have reproduced the numerical results of refs.0, = 0, and those 
of ref.||, = ~G V /V8, 5G 1 = 0. Nevertheless, the I (a, b) function in ref.Q was evaluated with the 
mass of the K° although the mass of the K + should have been used since it corresponds to a K + K~ 
loop. This kinematical source of isospin violation gives rise to non-negligible corrections due to the 
proximity of the K°K° threshold to the mass of the 0(1020). When this is taken into account, one has 
Br((j) -> -/K°K°) = 3.0 10" 8 instead of 5 10" 8 as given in ref.|?]] where the T-matrices of ref.|H]] were used. 
We now evaluate the previous branching ratio with our present formalism and with the values of the SG 
and Ci as given in table ||. For the strong amplitudes of ref. we obtain Br((ft — * jK K°) =3.7 10 -8 and 
with the T-matrices of ref. [34] one has Br(<p — > jK°K°) = 6.43 10 -9 . The tree level contact interaction 
(j)jK°K° and its iteration through final state interactions interfere destructively and tend to cancel each 
other or even they give rise to a negative interference with the pure K + K~ kaon loop contribution of 
ref.0. 



5 Conclusions 

In this article we have shown how the experimental results for the decay widths T((p — * 7/0 (980)) 
and T(0 — ► 7ao(980)) from refs.|l], ||, ||] can be described without abandoning isospin symmetry in 
the calculation of the S-wave strong T-matrix elements [31, 34]. Nevertheless, in our final results we 
have calculated the integral I(a, b) in terms of the K + mass instead of the average isospin mass. This 
kinematical isospin violating fact amounts to effects of around a 10% in the rate <j) — * 7/0 (980) and 
around a 20% in the width to 700 (980). The jK°K° branching ratio is much more sensitive to these 
kinematical effects due to the so much reduced available phase space. The same situation would have 
arisen in the decays (ft — * 7-R as well if we had used the standard two body decay formula eq. ( |3.10| ) , with 
well defined masses itir, instead of having taken care of the finite width effects of the /o(980) and ao(980) 
resonances. These results are quite opposite to what has been claimed in ref.|23j regarding the necessity 
of including large isospin violating effects in the couplings of the /o(980) and ao(980) scalar resonances 
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to the K + K~ and K°K° channels due to the proximity of the KK threshold to the nominal masses of 
the /o(980) and ao(980) resonances. We have also stressed that one should abandon the standard two 
body decay formula, with well defined masses for the /o(980) and ao(980) resonances, in the calculation 
of the rates <j> — > 7/0 (980) and <j) ~~ > 7 a o (980). Instead, finite energy distributions, /r(Q ), have to be 
considered from the very beginning because of the dramatic changes in the resulting decay widths under 
small changes of the /o(980) and ao(980) masses as compared to their widths or as compared to the 
difference between the pole masses and the energy of the peaks in the S-wave 1 = and 1 T-matrices. 

It is also worth remarking that the formula derived in sec.|| for calculating the decay widths <j) — > jR, 
eq.dOf), together with eq. ( 3.171 ), can be also applied in the experimental analyses of the rates (f> — > jR 



as an alternative to the rightful one followed in refs.[|l], ||, |], ||, || from ref.[]T2[. In this way, one does not 
use, even as an intermediate step, the bare theoretical concept T(cp — ► 7^; \J Q 2 ) and incorporates the 
contribution to W 7 (Q 2 ) proportional to (j in eq.( |2.12 ). 



In order to describe the data we have included, beyond the K + K loop model of ref.[12|, self-gauge 
invariant vertices Vi(<jyyKK) with the KK pair in the isospin channel I. As a result of the fit to 
the experimental data of refs.|^, |f| [|, |6|, a non-vanishing V(4r/K°K°) local term has emerged while 
the corresponding V{(jyyK + K~) tend to vanish. It is also shown that the former plays an important 
role in order to reproduce the experimental data. We have also included this contribution to the <j) — ► 
jK°K° decay although its effects do not spoil here the conclusions of ref.0], where only the K + K~ loop 
contribution is included, that this branching ratio is negligible small and does not offer any significant 
background to study CP violation in a ^-factory like DA<£NE. 

It has been repeatedly stated that the study of the </>(1020) radiative decays to jR constitutes an 
important test to unveil the nature of the /o(980) and oo(980) resonances by comparing the resulting 
experimental data with the models and approaches present in the literature. Indeed, we see from eq.flglg) 
that the study of these decays constitutes an alternative source of experimental information on the /o(980) 
resonance since it is sensitive to Imagft^] which is not directly measured in tttt scattering data. For the 
ao(980) resonance the experimental data is much more scarce than for the /o(980) and hence having new 
precise data, as that of ref.Q, is of foremost importance. Our simultaneous study of the new and accurate 
data of H on (ft — » 77r°r/ together with that of (f) —* jtt°tt° 0, |H| has given a coherent reproduction of 
these data in terms of only two free parameters. This constitutes a step forward in the experimental 
verification of the strong T-matrices of refs.[34, ^l[ already successfully tested in refs.[31, 32, 34|. 
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